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1 The Heisenberg spin chain is considered in 4 model approximation. Quantum cor- 

rections to classical solutions of the one-dimensional cj) A model within the correspondent 
physics are evaluated with account of rest d — 1 dimensions of a d-dimensional theory. A 
quantization of the models is considered in terms of space-time functional integral. The 
generalized zeta-function formalism is used to renormalize and evaluate the functional 
integral and quantum corrections to energy in quasiclassical approximation. The results 
are applied to appropriate conditions of the spin chain models and its dynamics, which 
elementary solutions, energy and the quantum corrections are calculated. 



^ Introduction 

m 

Q\ ; There is a wide field of Heisenberg spin chain [TJ realizations, intensely studied as quantum in- 
tegrable system [2] and in the context of its static and dynamic properties in external magnetic 
field [21 [3]. The original approach of W. Heisenberg is based on a localized electrons as initial 
CN ■ approximation, valid for metals with weak conductance. The second approximation accounts 
quantum exchange (due to Pauli principle) between electrons in different places. Starting from 
Heitler-London formula for the exchange and Coulomb integrals, one arrives at a (Heisenberg) 
^ ■ Hamiltonian which, by construction, describes spin system field of a solid. Investigations of 
symmetry in the first paper of Heisenberg [1] gives a fundamental approach to magnetics classi- 
fication as well as to ferromagnetism phenomenon and, for example, its existence only in cubic 
crystals (eight neighbors necessity). All this allows to believe in further developmenmt as the 
whole model as its particular cases applications. 

Some results of the known model applications are directly related to experiment in thermo- 
dynamics context, e.g. in [I]. Easy-plane ferromagnetism and in-plane domain-wall form factor 
is theoretically studied in connection with neutron scattering effects in CsNiF^ crystal [5]. Let 
us stress its quantum origin, which apart from well-known linear quasiparticles, provides a way 
to account for nonlinear collective phenomena as kinks, solitons and cnoidal waves. Generally 
speaking we have important applications of soliton theory aspects, as, e.g. in [B|. Some papers 
reduce the theory to Sine-Gordon model [7j. Other option leads to non-integrable 4 model, 
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which solutions are very similar. We would note a growing interest to a secondary quantization 
of such nonlinear quasiparticle fields, that, for example, allows to obtain so-called quantum 
corrections to classical energy of the objects. 

There is a method that is convenient for evaluation of the corrections. It is Feynmann 
integral by trajectories [5] (continual integral). The quantum corrections is a topic of stable 
interest since the seminal paper of R.F. Dashen, B. Hasslacher and A. Neveu [9], see also L. 
D. Faddeev, L.A. Takhtajan and V. E. Korepin papers [10]. The functional integral method 
becomes a practical tool for evaluation of quasiclassical corrections to the action from the times 
of V.P. Maslov paper [TTj . 

One of principal results of the method is obtained in [12], where the general algorythm of 
corrections evaluation is elaborated for arbitrary background profiles, expressions for ground 
state energies were derived for a 3+1 dimensions theory with a potential dependant on a single 
variable. Generalization for the supersymmetric kink is given in [13]. Solutions for Sine-Gordon 
quasiperiodic potentials was given in [T3] , 

Developing these results to arbitrary dimensions, investigating kink models and periodic 
solutions we demonstrated details of the Feynmann integral construction and generalized zeta 
function evaluation as well as the renormalization realization [7JIII2- A general algebraic method 
of quantum corrections evaluation based on zeta-function [16] is used and the Green function 
for heat equation with an elliptic potential is constructed (see also [T7]). 

In this paper we continue our investigations of the problem in the spirit of [7] and fix our 
attention on Heisenberg chain model in an so-called (J) 4 or Landau- Ginsburg [IB] approximation 
(see also Gross- Pit aevski equation [TjJ]) approximation. The merit of the presented paper is the 
attempt to derive conditions of possible application to realistic conditions of a magnetic medium 
(Sec. 2). In the next section we describe general features of the Heisenberg spin chain model 
and its reduction in specific conditions of 4 . In Sec. 3 we reproduce formulas resulting form 
[7J for readers convenience, namely, the space-time consideration, close to original Feynmann 
papers [5] . As for terminology, we use the word renormalization [7J to exclude divergence terms 
while in some other papers the term regularisation is met. The final section is devoted to 
specific case of so-called zero "mass" m 2 = condition, that is specified by a distinguished 
value of magnetic field, in which the field configuration drastically changes. 

2 The (f) 4 model of Heisenberg spin chain 
2.1 General equation of motion 

According to 0] (with ~t n = {3*, S^) as unit vectors) the Heisenberg magnetic chain with 
anisotropy in the direction of the chain and external magnetic field perpendicular to the chain 
is described by a classical Hamiltonian 

H = - JJ2 ■ + D J2(Sn) 2 ~ 9VbB £ S£ (2.1) 

n n n 

with corresponding equation of motion (here in SI units) 

Hdtt n = ~3 n x (- J{S n+1 + + 2DS*J - gn B Bx), (2.2) 

where J, D are spin coupling constants and g is the effective electron g-factor, Hb is the Bohr 
magneton, x and z are versors and B is the magetic field. After taking the continuum limit 



2 



(with a as lattice constant) one obtains 

hd t S x = -Ja 2 {S y d 2 S z - S z d 2 z S y ) + 2DS V S Z 
hd t Sy = -Ja 2 {S z d 2 z S x - S x d 2 S z ) - 2DS X S Z - g^ B BS z . (2.3) 
hd t S z = -Ja 2 (S x d 2 Sy - syd 2 s x ) + g^ B BSy 

By substituting S = (cos 9 cos 0, sin 9 cos 0, sin 0) one can reduce equations of motion to 

hcos<pd t 9 = Ja 2 (d 2 (f) + sin cos 4>{d z 9) 2 ) — 2.Dcos0sin0 — gfi B B sin cos 9, . . 

hd t (j) = —Ja 2 (cos(pd 2 9 — 2smcj)d z 9d z (j)) + g/i B Bsm9. ^ ' ' 

2.2 Model 4 approximation 

Stationary points of the system for D < (easy axis anisotropy) are such pairs (9, 0) for which 
= 0, G { — arccos ( ) , , arccos ( ) } with = being unstable. For g^ B B close to 
—2D both stable points are around = 0. In such a situation it is valid to assume 0^0 (with 
3 as the highest considered term) and with 9 as the highest considered term. 



hdtO = Ja 2 d 2 z <P - 2D(0 - 2f!) - gilB B{<j> - % ), 
frd 4 = -Ja 2 d 2 9 + gfi B B9. 

If we additionally assume | Ja 2 d 2 9\ << \gfi B B9\, we obtain 

hd t 9 = Ja 2 d 2 z (j) - 2D(<f) - - gfi B B((f) - f ), 
hd t (j) = gn B B9, 

which leads to 

d 2 <P = Ja 2 d 2 (j) - (2D + gfi B B)<p + 8g+ f B V , 



= ^a t c 



(2.5) 



(2.6) 



(2.7) 



ft — 

9^ B B 

The result represents the 4 model with the energy density 

h 2 (d4>\ 2 J a fd<f>\ 2 2D + gf i B B ^ SD + gfisB . 
H= 2^B{m) + TUJ + 2a " 24a ■ (2 " 8) 

After rewriting the equations in dimensionless variables (z = az' , t = Tt' with T as the time 
scaling parameter in the Feynman integral as in [7]) we obtain 



r4^d 2 4 = Jd 2 ^ - {2D + 9llB B)<j> + MZd^^s 

TT _ h 2 (d<t>\ 2 , J(d<p\ 2 , 2D+g/i B B ±2 8D+ gf i B B i4 
n ~ 2g/i b BT 2 {df) 2 \dz'J 2 V 24 ^ 



(2.9) 



For simplicity of further calculations we will write 

V 2 6(2D+g^ B B) 

8D+ 9 ij, b B 
m 2 _ 2D+gfi B B 

2 _ JgitbBT 2 
c ~ h 2 



(2.10) 



where m and V are parameters of the potential and c is a dimensionless propagation speed. 
Energy density takes the form 
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2.3 Static kink of the 4 model 

For the above described system there exists a well known static kink solution 

(f) = Vtimh(mz') , (2.12) 

which in this case represents a crossection of a flat, uniform domain wall in direction of it's 
normal vector. Classical energy of the kink is given by integration of the energy density (12. 8 p 

E = llJmV 2 = Uy/=J(2D + gii B B)l 
12 2V2(8D + gfi B B) 

For a domain wall this represents the energy per single chain of atoms. It is of note, that kink 
solutions vanish, when g\iB^ reaches —2D. There are however static solutions, which are still 
present for gfisB > —2D and an example will be discussed in section OU). 



3 Quasiclassical Quantum Corrections 
3.1 Quantization Scheme 

For the purpose of this publication we will use a semi-classical quantization procedure explained 
in detail in [TJ. For a given classical system described by action integral S with a static solution 
4> we derive energy corrections by expanding the action in path integral formulation of the 
propagator 

((f)\e-^ TH \(f)} = f D<p(x,t)ei s M (3.14) 

in a Taylor series around the classical solution and cutting it at the first non-trivial term. Then 
the formal expression for quantum correction to energy is 

A£ = -Ahi(det[L]), (3.15) 

where L is the second derivative of the classical Lagrangian up to a multiplicative constant 
arising from the gaussian integrals during the derivation process (see [11] and [7]). We use 
zeta-function renormalization scheme to deal with emerging infinities 

AE = -— lim — / r s ~ l (g L (r,^ - g Lo (r,^ \^))d^dr, (3.16) 

where it covers all variables of the classical system, gi is the Green function of the heat equation 

(J- + Lj g L (r, ^ ) = 6(r)6tf - ^ ) (3.17) 

and Lo is an operator analogous to L with a constant potential (representation of vacuum). 
Additionally, mass scale is used to cut logarithmic divergence in all relevant parameters (see 
and [7J). Often the intermediate steps of (I3.16P are defined explicitly as 



7(r) 



J(g L (r,^,~^) - g Lo {T,~£ ,~£))d~£ (3.18) 
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and 

COO = fH) . 







It was shown in [T7], with important link to generalized zeta function theory, that if L can 
be written as a sum of operators acting on independent variables L = Yli-Li, heat equation 
Green function for L can be written as a product of Green functions for Lj. We are using this 
property to account for arbitrary number of spatial variables of the classical system. 

3.2 Quantum Corrections to 4 kinks 

We now procede to calculate quantum corrections for energy using the above described gener- 
alized zeta function renormalization scheme with following form of operators: 

Lx — A - Am 2 + 6m 2 sec/i 2 (mz')^ , (3.20) 

A d 2 

= (3 ' 21 > 

A - J^- < 3 ' 23 > 

where d is the total number of spatial dimensions, A^-i covers all spatial variables except for 
z', I is the range of all additional spatial dimensions in multiple of a (due to the same rescaling 
as for z) and r is the mass scale. Laplace transform of the Green function diagonal for L\ was 
derived by use of algorithm described in [TB] . For L 2 and L 3 spectra continuum approximation 
was taken. We obtained following corrections 

AEd=1 = W^{ 2 + ^~ 2 ln(2) " 3 ln (-^n 2 )) , (3-24) 

hcm 2 l ( 3 / 1 \\ . 

kE d=2 = — — - 3 + - arcsin [ — )) : (3.25) 



2T7T V 2 Vv/3 

Hcm 3 l 2 f 2 \ 

AE d=3 = ( ~ 6 H-Am 2 ) - 6 + -(-11 + 3^ + 6 M(2))j . (3.26) 

Mass scale r is chosen so that any logarithmic contributions will vanish (ln(— Am 2 ) = 0) for 
arbitrary T. Energy corrections exhibit a similar dependance on classic equation parameters as 
those of Sine-Gordon system [7]. Physical meaning of those parameters is however different. 
It is of note, that similar results were obtained by Konoplich in [TB]. Difference comes from 
accounting for {^) term in classical action. After inserting proper forms of c and m we obtain 



=1 " 72V 1 1 + 27! " H ) ) ' ( 7) 

Vg7^B(2D + gt i B B)l ( 3 . / 1 \\ 

=2 = 477^ V + 2 arcsm l7I J J ' (3 - 28) 
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AE d - 3 = V=«Smo + S^)iP f_ 38 , 2 A 

In this system the difference between one-dimensional model and one accounting for two addi- 
tional spatial dimension is especially visible near the qubB = —2D border case. Since AEd=i 

1 3 

is proportional to (2D + g^sB)^ instead of (2D + gixsBp, it would outweight the classical 
energy significantly. If we look at the ratio of corrections to classical energy for d = 3 

AE^ = «(8P + g , B B) ( 38 , 2 A 

44ir 2 ji V 9 3 k ') y 1 

we can see a particularly strong dependance on J parameter, which represents interaction 
strength between neighboring electrons. Calculated ratio will be the highest for gfisB — > —2D. 
In this limit we obtain 

Afe f-D\i 

— K -(~J ' (331) 
It is worth noticing, that due to the way Planck constant enters propagation speed c, energy 
correction are independant of it's value. Regardless of the number of spatials dimensions taken 
into account, energy corrections still vanish along with the classical field, when gusB = —2D 
(see the next section). 



4 The special case of 2D + gfisB = 

In this section we will discuss a static solution of 4 model, which does not vanish, when both 
stationary points of the potential coincide. Let us rewrite equation of motion and Hamiltonian 
in such case, when m 2 = 2D + g^sB = . 

id 2 t 4 = Jd 2 z ,<j) + D0 3 

H=i(^m 2 +m 2 -ii^ • (432) 

This system has an interesting traveling wave solution 



<P(b(z' ± vt 1 )) = b\ 2J ^* f^ snibiz! ± vt')- 1) (4.33) 
V —Dc 2 

with b as a free parameter real-valued amplitudes for v < c due to D < (sn denotes Jacobi 
SN elliptic function). For now we will focus on the v = case 



b\l^sn(bz'-i). (4.34) 

It is particularly interesting, since it's a non-trivial static solution in a system with a single 
stationary point = 0. Quantum corrections to the energy of this solution can be calculated 
through the same procedure as for a (p 4 kink with 

L 1 = A^ + &8n 2 (bJ)Y (4.35) 
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where A is the same as in (I3.23P and other Lj as before. Let us reproduce some formulas from 
|17j . Laplace transform of the Green function diagonal can be obtained in the same way as for 
d> 4 or Sine-Gordon kinks [15] 



G l{ p, Z ) = P 2 -^Ml-z) + 9^-2)z (436) 
2v/(36 2 + p){3b 2 - p)p{p 2 - 126 4 ) 

It is important to note, that for calculation convenience, we use rescaling r — > ^ for the Green 
function equation ( 13. 17ft . The polynomial in the denominator has five simple roots. We will 
set the vacuum counterpart to coincide with the highest one 

G (p,z) = — 7 J==. (4.37) 

We can now integrate the Green function diagonal over the period , where K is the complete 
elliptic integral of the first kind 

y(p)= I " (Gi(j9,cn 2 (fox)) -G (p,cn 2 (bx)))dx, (4.38) 
lip) 



o 



6b 4 K(i)+2p 2 K{i)+36b 4 (K(i)-E(i))-3b 2 p(E{i)-3K(i)) 

by/{3b 2 +p)(3b 2 -p)p(p 2 -12b 4 ) . . 

_ 2K(i) ■ l 4 - dy J 

by^2tiW^p 

To properly define the inverse Laplace transform, we need the 7 function to be smooth in 
a oi < 3?(p) < o p area with 01 < o p as constants [20]. Since the 7 function has five distinct 
singularities (all on the real axis), we have six potential inverse Laplace transform (disregarding 
the choice of the sign for the square roots), but only one of them fulfills the necessary condition 
V T<0 7(7") = arising from the definition (13. 1 Tj) - it occures, when c p — > 00. Thus the 7 function 
as defined in (13. IS)) will take form 

ro+ioo 

7 (r) = / <? AT i(p)dp, (4.40) 



o— 100 



where > 26 2 \/3. At this point, we add all relevant variables as for 4 case 7(7-) — > 7 (r) 72(7") 73 (r) 
The energy corrections per wave period are finally obtained by performing the Mellin transform 
and taking the derivative at s = 

fc o poo ro+ioo 

AE = -—Um 7r r-SW-ysCr) / e^(p)dpdr, (4.41) 



Plugging expressions from (14.391) and relevant forms of 72 and 73 (see [7]) yields 
AE = _ ^-^vOT^- _i n m JV-l-i 

f e fSJ f 6fe 4 X(i)+2p 2 K(j)+36b 4 (K(i)-g(0)-3b 2 p(g(0-3^(Q) _ 2AT(Q A rf , 

J °- io ° I 6 a/(36 2 +p) (36 2 -p)p(p 2 - 12b 4 ) by / 2b 2 V3-p J 



(4.42) 



There still remains the problem of finding the explicit analitic solution of the shown integrals. 
Applications of the distinguished condition of the elliptic field configuration may be interesting 
from the point of measurements realization. Such configuration may be more easily noticed 
(recognized in experiments). 
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5 Conclusion 



Energy corrections for easy axis domain walls would be particularly interesting in the case of 
thin ferromagnetic films, where they should dominate the classical energy. Formally, the results 
of the last section partially coincide with the known case (see e.g. [T7], where the so-called Nahm 
model of Yang- Mills theory is studied) but physical sense is different. 

The space-time consideration we develop in our publications has obvious intensions to in- 
clude mutisoliton configuration into quantum quasiclassical picture. 
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